Trapping molecular ions formed via photo-associative ionization of ultracold atoms 
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The formation of ^"Ca^ molecular ions is observed in a hybrid *°Ca magneto-optical and ion 
trap system. The molecular ion formation process is determined to be two-photon photo-associative 
ionization of ultracold ''"Ca atoms. A lower bound for the two-body, two-photon rate constant is 
found to be /3 > 2 ± 1 X 10~^^ cm"^ Hz. Ah initio molecular potential curves are calculated for 
the neutral Ca2 and ionic Ca^ molecules and used in a model that identifies the photo-associative 
ionization pathway. As this technique does not require a separate photo-association laser, it could 
find use as a simple, robust method for producing ultracold, state-selected molecular ions. 



The low-energy internal structure of a diatomic 
molecule, e.g. the electric dipole moment and vibrational, 
rotational, and fi-doublct levels, presents a host of op- 
portunities for advances in quantum simulation, preci- 
sion measurement, cold chemistry, and quantum infor- 
mation As such, the last decade has witnessed an 
enormous effort towards producing ultracold molecules 
in well-defined rovibrational states. It appears that this 
goal is in now within reach, as several groups have re- 
ported production of ultracold molecular samples in the 
rovibrational ground state 043 ; in the case of the KRb 
work the sample is near quantum degeneracy 0] . 
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To date, most work has focused on neutral molecules, 
however, a sample of ultracold molecular ions presents 
interesting possibilities. While the Coulomb interaction 
overwhelms the electric dipole interaction between po- 
lar molecular ions in most cases, many of the goals of 
cold molecule physics can be accomplished with molecu- 
lar ions, with the added benefit of a simple, reliable trap- 
ping. Indeed, a cold sarnple of molecular ions allows the 
study of cold chemistry [5[ , which not only has important 
implications for understanding the formation of interstel- 
lar clouds 3 1 but for investigation, and possible control, 
of reactive collisions in the quantum regime 0; identi- 
fication of carriers of the diffuse interstellar bands Q; 
precision measurement of molecular transitions, which 
can be used to sensitively measure parity violating ef- 
fects 0, as well as to constrain the possible variation 
of the fundamental constants [lo| : and the implemen- 
tation of a scalable quantum computation architecture 
llj. For these reasons, several groups have recently ini- 
tiated work to produce samples of ultracold molecular 
ions [iMi- In Refs. [i3-[ii, the molecular ion exter- 
nal degrees of freedom are sympathetically cooled with 
laser cooled atomic ions to ultracold temperatures. How- 
ever, the long range ion-ion interaction does not couple to 
the molecular ion internal degrees of freedom and there- 
fore the rovibrational temperature is unmodified by the 
sympathetic cooling process. To overcome this limita- 
tion, optical pumping schemes [13, [l^ and state selec- 



FIG. 1: MOTION trap: A ■*°Ca magneto-optical trap is con- 
structed inside a linear quadrupole ion trap. Below the trap 
are two ablation targets used to load calibration ions. The 
channel electron multiplier (GEM) is mounted above the trap. 



tivc photo-ionization of neutral molecules jl4| have been 
used to produce molecular ions in the lowest few rota- 
tional states. More generally, it should also be possible 
to use ultracold neutral atoms, in a magneto-optical trap 
(MOT) , to simultaneously sympathetically cool both the 
internal and external degrees of freedom of molecular ions 

In addition to these methods, it is possible to create 
state-selected ultracold molecular ions directly from their 
atomic constituents through photo-associative ionization 
(PAI), a two-body process where colliding neutral atoms 
arc photo-excited into either a bound molecular ion or 
auto-ionizing molecular potential. In this manuscript, 
we report the formation of '^^Csl^ molecular ions via PAI 
in a hybrid system composed of a '"'Ca MOT constructed 
inside of a radio-frequency linear ion trap (MOTION 
trap) . We quantify the production of the dimer ion using 
standard ion trap mass spectrometry techniques 
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FIG. 2: Results of a excitation spectrum measurement per- 
formed at trapping parameters where the stabihty parameter 
(q) for m = 80 amu is 0.5, thus unstable for m = 40 amu. 
The observed loss near 70 kHz is within 10% of the predicted 
secular frequency for m = 80 amu. The split peak feature 
is likely due to trap imperfections (asymmetric endcap place- 
ment, machining, imperfect quadrupole ratio, etc). 



available for use in the MOTION trap system. In the 
remainder of this work we briefly describe the apparatus 
used and present measurements characterizing the molec- 
ular ion production rate. We conclude with results from 
ab initio calculations of the ^"^Ca^ molecular ion struc- 
ture and a description of the most-likely mechanism of 
molecular ion production. 

The MOTION trap system used in this work is that 
proposed in Ref. 15| and is shown in Fig. [TJ As afore- 
mentioned, it is a hybrid system with a ^°Ca MOT con- 



structed inside an linear quadrupole ion trap (LQT), de- 
signed to explore collisions between ultracold atoms and 
trapped molecular ions. Briefly, cold ^''Ca atoms are 
collected in a MOT by decelerating a beam of atoms 
produced by a Ca getter. The A = 422.6 nm laser sys- 
tem , used to produce the MOT cooling and decelera- 
tion beams, is referenced to the 4^Pi <— 4^So transition 
(r = 277 X 34.6 MHz) using a saturated absorption signal 
from a Ca vapor cell. The trapping (deceleration) laser 
beams typically have detuning and total intensity of — F 
(— 3F) and 23 (31) mW/cm^, respectively. The 671.7 nm 
5^Pi <r- 3^D2 Ca repump transition is driven with a laser 
system referenced to the R69(2-5) molecular Iodine tran- 
sition. The typical MOT atom number, density and tem- 
perature are measured by absorption and fluorescence 
imaging and found to be iVca = 1-0 ± 0.1 x 10'' atoms, 
pCa = 3.9 ± 0.5 X 10^" cm-3, and Tca = 4 ± 1 mK, re- 
spectively. The LQT, which has a field radius of ro = 
11.2 mm, is designed with three axial trapping regions 
to facilitate the shuttling of ions into and out of the cen- 
ter region for interaction with the MOT. The number of 
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FIG. 3: Measurement of the ion number as a function of load- 
ing time is shown in panel (a). In panel (b), trap lifetime is 
measured by loading the trap to saturation and stopping pro- 
duction for a fixed amount of time before detection. Charac- 
teristic times are found by exponential fits (dashed lines). 



trapped ions is detected by grounding the trap electrodes, 
allowing the ions to escape and be detected by a chan- 
nel electron multiplier (CEM) ion detector located near 
the trap center. The absolute ion number measured by 
the CEM is calibrated by comparing CEM detection to 
laser-induced fluorescence rates with ^^^Yb+ ions loaded 
via laser ablation (see Fig. [IJ. Only the central trapping 
region is operated in this experiment since the MOT acts 
as the source of all ions observed in this work. 

When the MOTION system is operated without de- 
liberately loading ions into the LQT, an accumulation 
of ions is observed due to the presence of the MOT. In- 
terestingly, though there exist several pathways to the 
formation of "'^"Ca"'' - e.g. photo-ionization of an atom 
in the 5^Pi state, populated via the repumping chan- 
nel, by a 423 nm photon - an ion signal is still observed 
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FIG. 4: Ion number measured as a function of the MOT den- 
sity determined by absorption imaging. A quadratic trend 
(with a best fit shown to guide the eye) is observed until trap 
saturation occurs at higher densities (shown in the points with 
open circles). 



even when the LQT is operated such that the charge-to- 
mass ratio of ^'^Ca"'" is not trapped. The mass of this 
anomalous ion species is measured by the LQT stability 
parameters, a comparison of the relative CEM time-of- 
flight to that of ions with known mass {e.g. *°Ca+) and, 
ultimately, by resonant excitation of the secular trap fre- 
quency as shown in Fig. [21 The excitation measure- 
ment is performed by applying a small (relative to the 
trapping potentials) oscillating voltage to one of the trap 
electrodes, inducing substantial heating that results in 
trap loss when the excitation is resonant with the mass- 
specific secular frequency of the ion motion. The excita- 
tion spectrum is taken by allowing the LQT to load with 
ions for 1.9 s at each excitation frequency before detec- 
tion with the CEM. The ratio of detected ion number 
with excitation to detected ion number without excita- 
tion is plotted in Fig. [2] The trap loss resonance (and 
its second harmonic) is observed as the reduction in ion 
signal near 70 ± 2 kHz (140 ± 5 kHz), in good agreement 
with the predicted radial secular frequency of 76 kHz 
for 80 amu, confirming the production of ''°CaJ molecu- 
lar ions. The splitting of the fundamental peak is likely 
due to trap imperfections such as the asymmetric axial 
electrode placement, asymmetric voltages, and imperfect 
machining. 

The LQT loading dynamics of the dimer ion arc shown 
in Fig. O For the data in Fig. I^l^a), the LQT is first 
loaded with dimer ions from the MOT for a designated 
amount of time after which the ion number is measured 
by the CEM. As seen in this panel, the characteristic 
loading time is found to be ^ 1 s and typically leads to 
an equilibrium ion number of ^ 3000 ions. To measure 
the competing effects which lead to this equilibrium, we 



measure the trap lifetime as shown in Fig. [Sfb). For this 
measurement, we allow the trap to load to saturation 
and then switch off the MOT at t = s, stopping molec- 
ular ion production. After a designated time, we detect 
the ion number via the CEM. In this way, we find the 
molecular ion lifetime in the LQT to be ~ 90 s. Since in- 
elastic processes between the ions and neutral atoms [l3| 
do not change the number of dimer ions in the trap, the 
loading time is not affected by the presence of the MOT. 
Therefore, the large disparity between the molecular ion 
lifetime and its characteristic loading time shows that 
the equilibrium condition is not set by competition be- 
tween the natural trap loss mechanisms and the dimer 
ion production rate. Instead, equilibrium is determined 
by a maximum ion density that sets a trap capacity at 
which the Coulomb repulsion between the ions already 
in the trap and the freshly formed ions is large enough to 
overwhelm the trapping potential [l7| . Further evidence 
to support this conclusion comes from measurements at 
different MOT parameters with different characteristic 
loading times which reach the same approximate equi- 
librium ion number, as well as measurements counting 
"'"Ca+ and ^^"'Yb+ loaded by photo-ionization and laser 
ablation, respectively. 

To confirm that the dimer ion is produced by collisions 
within the MOT and not by other means, such as colli- 
sions of MOT atoms with hot Ca from the getter, we 
measure the production of molecular ions as a function 
of MOT density, as shown in Fig. |4l For the data in 
this figure we detect the dimer ion number after 1.9 s 
of loading at various densities. We vary the density by 
controlling the intensity of the 672 nm repump laser; the 
density at each point is determined by absorption imag- 
ing. Quadratic dependence is found for sufficiently small 
densities, but trap saturation, due to the space charge 
effect, leads to constant ion number at higher densities. 
The dependence of the molecular ion production rate at 
smaller densities is indicative of a two-body process as 
would be expected for PAI between MOT atoms. 

To determine the number of photons required for the 
PAI, we conduct a similar experiment as a function of 
overall 423 nm laser light power with results shown in 
Fig. [5l Since the position and density of the MOT are 
sensitive to the exact power balance between the decel- 
eration and trapping beams, there is higher systematic 
uncertainty in this measurement. However, we observe 
a trend that is consistent with quadractic dependence, 
drawn as a dashed line to guide the eye in Fig. \E\ Thus, 
we conclude that the ^''Ca^ molecular ions are most likely 
produced by a two-body, two-photon production mecha- 
nism as detailed in what follows. Because space charge 
effects dominate the LQT loading dynamics, the molec- 
ular ion production rate cannot simply be analyzed in 
terms of the competing rate equation model typically em- 
ployed in trap loading dynamics. Therefore, in order to 
estimate this rate, we postulate that the number of Ca^ 
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FIG. 5: Normalized production rate as a function of overall 
423 nm laser power incident upon the MOT. A quadratic fit 
is shown as the dashed line. 
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FIG. 6: Molecular potentials of *°Ca2 (black curves) and 
^°Ca^ (red curve). Two vertical arrows indicate the energy 
carried in the 423 nm MOT photons. We observe that the 
ground state potential of the molecular ion is in the same 
energy region as the potentials dissociating to two '^P *''Ca 
atoms. 



molecular ions, TV, in the ion trap as a function of time 
is given as: 
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where the average production rate is characterized by 
/3, the two-body intensity dependent rate coefhcient; f] 
is the average of a spatially varying capture efhciency 
which quantifies the probability of successfully trapping 
an ion formed at a particular location; and V are 
the average squared MOT density and volume, respec- 
tively. The last factor is an ad hoc model to take into 
account the trap saturation at the observed maximal ion 
number No ^ 3000 ± 300. This model allows us to use 
full loading curve measurements to set a lower bound on 
the production rate in a way that is equivalent to per- 
forming a model independent linear fit to only the very 
early {N ^ ^ 0) portion of the loading curve, where 
space charge effects should be minimal. The rate coef- 
ficient from the fit, assuming perfect trap capture effi- 
ciency [fj — > 1), to the loading curve shown in Fig. [3] is 
found to be /3 = 2 ± 1 X 10~^^ cm^ Hz. This rate coeffi- 
cient should only be interpreted as a lower bound to the 
true rate constant since the capture efficiency is taken to 
be its maximal value. Interestingly, this rate coefficient 
is comparable to observations in alkali systems [l^, 
when scaled to the relatively low intensities found in this 
experiment. This bound is independent of any assump- 
tions about the specific mechanism leading to the PAI 
event. 

In order to determine the specific PAI pathway for 
molecular ion production, molecular potential curves for 
the neutral ^°Ca2 and ionic ""^CaJ molecules are pre- 
sented in Fig. m The potentials of the XiE+ and B^Y.+ 
states of ^°Ca2 are constructed from spectroscopic data 
in Refs. [U 22 1. Other potentials are determined us- 



ing non-relativistic configuration-interaction molecular- 
orbital restricted-active-space (MOL-RAS-CI) calcula- 
tions. From these molecular potentials it is clear that 
there are several energetically allowed pathways for asso- 
ciative ionization that must be considered: 



Ca + Ca+ 
Ca(P) + Ca(S) +7 
Ca(P) -I- Ca(P) 
Ca(D) + Ca(S) +7 
Ca(S) + Ca(S) + 27 



Ca.2 
CaJ 
Ca^ 
CaJ 
Can 



-7 
e~ 



(1) 
(2) 
(3) 
(4) 
(5) 



The first pathway ([ij is the radiative association of the 
atomic ion with its neutral partner. To investigate this 
reaction, we trap large amounts of ^"Ca"'' and monitor 
the production of ^°CaJ in the presence of the MOT. 
We see no evidence of this production mechanism and 
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conclude that radiative association cannot be the domi- 
nant pathway. 

The next two pathways are single step processes in 
which one of the colhsion partners resides in an excited 
state before the colhsion reaches the Condon point, where 
a free-to-bound transition is possible. The excited colli- 
sion pathways (|2l3p are ruled out because the duration 
of the cold collision is longer than the excited 4^Pi state 
lifetime. As discussed in [l9[, excited atoms are unlikely 
to be in close enough proximity to undergo association 
since the 423 nm laser driving the transition is far de- 
tuned due to the inter-atomic interaction even at dis- 
tances large compared to the Condon length. 

The other excited collision pathway ([4]) cannot be ruled 
out due to collision duration argument [l^ because the 
3^D2 Ca state is relatively long lived, with a radiative 
lifetime ~ 1 ms, and could survive for the duration of a 
collision event. To quantify this production mechanism, 
we must consider the dependence of the dimer ion pro- 
duction on the partial density of 3^D2 Ca in the MOT. In 
steady state, the density of 3^D2 state Ca can be deter- 
mined as PD = (TpD/OCaPee) / (Frad + Tke) where pCaPee 

is the density of atoms in the 4^Pi state, FpD is the ra- 
diative decay rate for the 4^Pi 3^D2 transition, Fj-ad 
is the total radiative loss rate from the 3^D2 state, and 
Fke is the ballistic escape rate from the MOT volume. 
The ballistic escape rate is determined to be 1.4 ±0.2 kHz 
by fitting to the cloud volume evolution measured during 
free expansion of the MOT atoms. The total radiative 
loss rate is determined by solving a rate equation model 
for the effect of the repump laser on the 3^D2 state pop- 
ulation that includes the 5^Pi and 4^ So states as well as 
the losses into the 3'^P2.i states. Using this system of rate 
equations, the partial density of 3^D2 ''^Ca is determined 
for each data point in Fig. |4l The result, after dividing 
out the overall quadratic density dependence, is plotted 
in Fig. [71 Despite the fact the density of 3^D2 Ca varies 
by over factor of two, the rate constant is unchanged. 
We note the two points in open circles correspond to the 
points in Fig. |4]where the ion number is space charge lim- 
ited, and should be ignored. Therefore it appears that 
the 3^D2 is unimportant for the observed PAI and thus 
pathway ^ is ruled out. 

Therefore, we conclude that the two-photon PAI pro- 
cess (O is responsible for the observed production of 
^'^Ca^ molecular ions. This formation proceeds by the 
excitation scheme sketched in Fig. |6l where colliding 
neutral atoms are photo-associated into an excited neu- 
tral molecule state, which is subsequently excited to the 
ground molecular ion state or a further excited, auto- 
ionizing neutral molecule state. We note that this path- 
way agrees well with the conclusions of Ref. 
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In summary, we have observed the formation of ^''CaJ 
molecular ions and found a lower limit of its produc- 
tion rate constant at ^ > 2 ± 1 x 10"^'"' cm^^ Hz. From 
the dependences of the production rate on MOT pa- 
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FIG. 7: Ion number normalized by the square of the MOT 
density as a function of 3^D2 state density. Note the open 
circle points correspond to those in Fig. |4] at which trap 
capacity is reached. 



rameters and ab initio calculations, we have determined 
the production mechanism is most likely two-step, two- 
photon photo-associative ionization out of the ground 
state molecular potential. We produce approximately 
3000 molecular ions, limited only by the space charge ca- 
pacity of the ion trap. As this technique does not require 
a separate photo-association laser, it could find use as 
a simple, robust method for producing ultracold, state- 
selected molecular ions. 
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